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1. Introduction 

The discovery of a new particle with a mass of about 125 GeV at the LHC [1,2] may be 
explained as discovery of the Standard Model (SM) Higgs boson. But it is also possible that the 
lightest MSSM Higgs boson h° has been found. The identification of the new particle as SM 
or MSSM Higgs boson requires precise theoretical predictions for its production and decay. A 
lot of work has already gone into the more and more improved predictions of production cross 
sections and decay rates, which goes hand in hand with including effects arising at higher orders in 
perturbation theory. A short review of the current status concerning the production cross sections 
of the Higgs boson has been given in the talk connected to this proceedings contribution. It was 
mainly based on the information given in the comprehensive reviews [3] and [4]. Such a review will 
be skipped here in order to be able to focus on the new contributions presented at the conference, 
which have been published in [5,6]. 

In detail this concerns the calculations of the total production cross section of h° via gluon 
fusion at next-to-next-to-leading order (NNLO) within the effective theory framework. In more 
detail this calculation is enabled by the determination of the effective Higgs gluon coupling C\ up 
to three loops in the strong coupling constant a s in the MSSM. 

In the SM it is well known that radiative Quantum Chromodynamics (QCD) corrections to the 
loop induced gluon fusion process do play an important role. The next-to-leading order (NLO) 
QCD corrections can increase the leading order (LO) cross section up to +100% [7, 8, 9, 10] 
including a large uncertainty due to the variation of the cross section in dependence of the renor- 
malization and factorization scale [l. The corrections appearing at NNLO further increases the 
cross section up to +30% [11, 12, 13] compared to the NLO predictions. They have first been 
calculated using the effective field theory approach, where the mass of the top quark is infinitely 
heavy. It turned out that after including these NNLO corrections the scale variation is strongly re- 
duced and one can obtain a trustable prediction using perturbation theory. By taking finite top mass 
effects at the NNLO level into account, it has been checked, that the used effective field theory 
(EFT) leads to reliable predictions [14, 15, 16, 17, 18, 19] as long as the Higgs mass stays below 
about 200 GeV. 

The calculations of the production cross sections of h° in the MSSM at NLO has been accom- 
plished by many groups using different methods including bottom induced contributions, which 
can be large compared to the SM ones [10] due to a tan/3 enhanced bottom Yukawa coupling [20, 
21,22, 23,24, 25,26, 27, 28]. 

Before the publication of Refs. [5, 6] the NNLO corrections had to be approximated by the 
corresponding SM corrections [21, 29], because the three loop part of the effective Higgs coupling 
C\ was not known in the MSSM. In a first step this missing piece has been calculated with degen- 
erated masses for supersymmetric particles (sparticles) in Ref. [5]. Ref. [6] extends this calculation 
to apply to more general supersymmetric mass spectra. 

The layout of this proceedings contribution will be as follows: In chapter two we shortly 
review the calculation of the effective Higgs gluon coupling at three loop order. In chapter three we 
show the connection between C\ and the production cross section of h°, which will be evaluated 
numerically in chapter four in certain scenarios at NNLO. The last chapter contains summary and 
conclusions. 
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2. Effective Higgs gluon coupling 



The effective Higgs gluon coupling arises when one integrates out all particles heavier than 
the Higgs boson H/h° (SM/MSSM) from the SM/MSSM Lagrange density. Thus in the SM one 
has to integrate out the top quark only. In the MSSM in addition all supersymmetric particles have 
to be integrated out. The only relevant remnant of those particles in the Lagrange density is the first 
term on the r.h.s. of the following equation: 



F.eff 



QCD- 



(2.1) 



where the Operator &\ is built up by the gluonic field strength tensor G°' /JV as follows: 

_n 1 



0,/xv 



(2.2) 



and -^qq D is the well-known five flavour QCD Lagrange density. To ensure the agreement of 
the predictions of the effective Lagrange density J2?y, e ff and the the supersymmetric QCD (SQCD) 
Lagrange density in the limit where the mass M/, of the Higgs boson is much smaller than the mass 
of the top quark and all supersymmetric particles, one has to determine the effective coupling Cj 
by comparing equivalent Green's functions of both Lagrange densities incorporating a single Higgs 
boson and gluons. As simplest Green's function one can consider the transition of two gluons to 
the Higgs boson like depicted for the SM case in Fig. l(a-b). The comparison of the two Green's 
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Figure 1: (a-b) LO Feynman diagrams for the process gg —> H in the SM. (a) In the full theory the Higgs 
boson couples via a top quark loop to the Higgs boson, (b) In the EFT a direct coupling of gluons and the 
Higgs boson emerges, (c) Particles appearing in the loops during the calculation of C\ in the MSSM. 



functions is performed in the on-shell (OS) production kinematics. That means (p\ + p-i) 2 
holds for the four momenta p\ and p% of the two gluons with p\ = p\ = 0. In order to determine C\ 
in the SM at LO one has to calculate the Feynman diagram in Fig. 1(a) (and the one with reversed 
fermion lines) in the Limit Me <C m t , where m t is the top quark mass. C\ is directly given as 
the first term in the corresponding expansion for a small Higgs masses and is proportional to a s 
at LO. For a more detailed description how to extract C\ the reader may be referred to Ref. [30]. 
In the MSSM case additional diagrams including sparticles have to be considered. A list of all 
particles appearing inside loops can be found in Fig. 1(c). For a determination of C\ at NLO/NNLO 
all relevant two/three loop diagrams contributing to this process have to be calculated. Example 
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diagrams for the MSSM case are shown in Fig. 2. Due to the expansion in the hierarchy M# <C m t 
the problem is reduced to the calculation of two/three loop tadpole integrals with more than one 
scale. At two loops the analytic results of the integrals for arbitrary masses are known [31] and C\ 



can be calculated exactly to this order. At three loops this is not the case and one has to reduce 
the number of scales in one diagram to one. This can be achieved assuming nearly degenerated 
and/or very different (hierarchic) masses. In the first case one can Taylor expand in small mass 
differences in the second case one can apply an asymptotic expansion. Combinations of both 
methods are possible, too. Once the number of scales has been reduced to one, the appearing 
integrals can be solved completely automatically with the FORM program MAT AD [32]. Because of 
the large number of diagrams, which have to be calculated, an automatic setup is mandatory. For 
the generation of Feynman diagrams we have used QGRAF [33], followed by q2e and exp [34], 
which provide an automatic application of asymptotic expansions. 

Besides the problem of multiple scales in MSSM diagrams, additional subtleties arise, when 
one uses Dimensional Reduction (DRED) as regulator, in order not to spoil supersymmetry (SUSY): 

Because the EFT cross section [35, 11, 36] was calculated using Dimensional Regularization 
(DREG), one has to ensure, that the epsilon scalars (which are the 2£-dimensional components 
of gluon in DRED) are integrated out, so that the EFT remains regulated in DREG. This can be 
achieved by giving the epsilon scalars a mass which is formally much larger than the Higgs boson 
mass. Due to the fact that SUSY is softly broken in the MSSM, the epsilon scalars do acquire a 
non vanishing mass via radiative corrections anyway. This mass has to be understood as evanescent 
parameter in the Lagrange density. 

In addition to that an evanescent coupling Ag of the epsilon scalars to the Higgs boson does 
appear via radiative corrections. This coupling has to be properly renormalized 1 at the two-loop 
level when one considers C\ at three loops. 

Apart from these subtleties one has to renormalize all particle masses for the three-loop cal- 
culation of C\ at the two-loop level, which are well known in the DR scheme [37]. Except for the 
gluino and epsilon scalar mass, where only the one loop renormalization constants are required. 

'for more details the reader shall be referred to Refs. [5, 6] 
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Figure 2: Feynman diagrams contributing to C\ in the MSSM. 
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To arrive at a finite result for C\ one requires the MS renormalization constant of the effective 
Operator Z\ \ , the decoupling constant of the gluon field £3 and the decoupling constants of GC S 
at the two loop level. 

Up to now C\ has been calculated at the NNLO assuming the following hierarchies: 
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Where the mass of the epsilon scalar was treated like being much smaller than the masses appearing 
above. The renormalization was done in the OS scheme, so that in the final expression of C\ the 
limit M e — > can be taken. 

To ensure that the described procedure consistently leads to C\ , it has been applied in the SM 
case. That means instead of using the common DREG regulator, DRED has been used like written 
above. Although the application of DRED to non supersymmetric theories leads to additional 
evanescent couplings, the well known SM (QCD) result was obtained when C\ was expressed in 
terms of the five flavour QCD coupling [5]. 

Another check for the calculation of C\ was enabled by the three loop calculation of £« t and its 
connection to C\ via the Low Energy Theorem (LET) [38]. It turned out that with the renormalized 
version of the LET one is able to calculate C\ from without even knowing the renormalization 
constant 5 A 2 of the evanescent coupling of the epsilon scalars to the Higgs boson. However, both 
ways lead, of course, to the same result for C\ . 

Explicit results for C\ can be found in Ref. [6] and Ref. [39]. 



3. Calculating Cross Sections 

The calculation of the production cross section of a light Higgs boson using the effective 
Lagrangian J2?y. e ff could be done [11] completely independent of the calculation of C\, because the 
effective coupling factorizes in the formula for the total cross section: 

a t ^) = a Q {^ ) ) 2 C\{^'^). (3.1) 
c i 

Where C\ = — ^~Y^=o (^r) c i"' an( * ^'0*) conta i ns the sum over all channels which in turn 
are obtained by convolutions of the partonic cross sections with the particle distribution functions 



(PDFs). Since Go contains the full Higgs mass dependence 2 at LO, Cbai 5 ^' 2 



gives the 

LO 

exact LO cross section. Thus in Eq. (3.1) higher order corrections obtained in the EFT are added 
relative to the factorized, exact leading order cross section, which improves the EFT prediction. 

A further improvement of this prediction can be achieved via a separation of hard \ih and soft 
scale \i s : It is known that only the combination C\(fl)ff\(n) is independent of the scale pL (up to 
neglected higher order corrections). Introducing the factor 3 B^ 5 \ which depends on the coefficient 



2 For the explicit formula in the MSSM see for e.g. Eq. (29) in Ref. [6]. 
3 The explicit definition is given in Eq. (32) Ref. [6]. 
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of the beta function and a. s enables the definition of a separately scale invariant matching co- 
efficient C g and effective operator & g , which do only depend on pL due to neglected higher order 
corrections. Thus one can choose a hard scale \ih in the matching coefficient and a soft scale \k s 
in the matrix elements of the effective Operator in order to avoid large logarithms. In that case the 
cross section formula becomes: 

Gt(Hs,Hh) = o C 2 g (n h )L(n s ) . (3.2) 

Where C g = J^Q and I = (fi^) 2 !'. Note that in the limit \i h = \i s the B-Factor cancels 

and one obtains Eq. (3.1). Eq. (3.2) holds for SM and MSSM predictions, as long as one uses the 
corresponding expressions for Ob and C g . However there is only one expression for L(n s •). Thus, as 
soon as the matching coefficient in the MSSM cj>QCD is known, it is straightforward to determine 
the cross section cr f S< ^ CD (jU,s, /i/,) in the MSSM using L(/J, s ) extracted from the SM cross section 
o® CD (jJ.s,IAh)- The predictions in the next section are based on the predictions of the SM cross 
section (jU 5 , jU^) obtained by the program Xsection, which was used to produce the SM 
predictions ofRef. [19]. 



4. Numerical Results 

In Fig. 3(a) the NNLO EFT cross section a ( S( ^ CD (^ s ,^/ ! ) for a hadronic center of mass energy 
of 8 TeV at the LHC is shown in pico barn assuming the wi™ ax scenario [40]. For the calculation 
the MSTW2008 PDFs set [41] and the result for cf QCD obtained in the hierarchy (hi) has been 
used. The Higgs boson mass M\ r between the red dotted (lower) line and red dashed (upper) line 
varies between 123 and 129 GeV. The strong decrease of the cross section for small values of Ma 
and large values of tan jS results mainly from a decrease of the top quark Higgs boson coupling due 
to the factor <=c which is also present in the corresponding top squark Higgs coupling. 
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The ratio of the SM EFT cross section and the MSSM one is shown in Fig. 3(b). Here the 
SM cross section depends only via the mass of the Higgs boson on the parameters Ma and tan /3 . 
That means for each point both cross section have been calculated using the mass of the lightest 
MSSM Higgs boson at this point. The gray shaded areas are already excluded by LEP [42] and 
CMS [43]. Only the parameter sets corresponding to the unshaded area between the two red curves 
are in agreement with current experiments. From the figure one can clearly see, that for growing 
values of Ma the cross section becomes more and more SM like. For small values of tan j3 and Ma 
a^ CD (iJ. s ,iJ.h) can be much smaller than a^ CD (iJ. s ,fXh)- 

At this point it is very important to note that in the used EFT approach no finite bottom mass 
effects have been considered. That means the bottom quark is treated as massless particle and 
does not couple to the Higgs boson. However for large tan /3 values (and small Ma values) it is 
know that the bottom Yukawa coupling to the Higgs boson can be enhanced. Thus one can expect 
large bottom quark induced contributions in this parameter region. In fact this is what other groups 
confirmed at NLO [24, 25, 22] and in order not to miss those contributions they have been included 
in Fig. 4(a). In addition the contributions of the known electroweak NLO corrections have been 

SOCD 

included. More details about how to combine a, with the mentioned contributions to obtain 
a SQCD can be found m Rfif [6] 

We want to note that the usual estimation of the NNLO part of c^ CD by the one of C^ CD 
results in a good approximation of the cross section in the m™ ax scenario. Where the error is clearly 
below the one percent level. 

In Fig. 4 a t (/J.) is shown in dependence of the common scale /I = pL s = \ih in the modified ra™ ax 
scenario, which was introduced in Ref. [6]. The dotted/solid curves represent the SM/MSSM result, 
where from bottom to top the blue/green/red lines show the LO/NLO/NNLO predictions. One can 
see a reduction of scale dependence when one increases the order of perturbation theory. The 
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observation that at NNLO there is a clearly visible difference between the QCD and SQCD curve 
makes clear, that the above statement concerning the estimation of the NNLO part of c^ CD given 
for m™ ax scenario does not hold for this scenario. In fact the scenario was constructed in such a 
way, that besides providing a Higgs mass in the current LHC range, we have a light top squark mass 
of about 370 GeV, whereas all other sparticles have masses around 1 TeV. This enables diagrams 
including the light top squark t\ to give sizable contributions to cf QCD , leading to a strong deviation 
from the NNLO contribution in the SM. 

5. Summary 

In this proceedings contribution a short introduction concerning the Higgs boson production 
via gluon fusion was given. Furthermore the calculation of the NNLO contributions to the effective 
Higgs gluon coupling cj !QCD in the MSSM was presented in a compact form. Selected numerical 
results implementing this contributions were shown. 

In the m™ ax scenario the SM NNLO part of cf CD can be used to estimate the NNLO part of 
£,sqcd However, it has been shown that in the presents of light sparticles like light top squarks, 
this estimation can fail. 
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